Exact solutions for equilibrium configurations of charged conducting liquid jets 
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A wide class of exact solutions is obtained for the problem of finding the equilibrium configura- 
tions of charged jets of a conducting liquid; these configurations correspond to the finite-amplitude 
azimuthal deformations of the surface of a round jet. A critical value of the linear electric charge 
density is determined, for which the jet surface becomes self-intersecting, and the jet splits into two. 
It exceeds the density value required for the excitation of the linear azimuthal instability of the 
round jet. Hence, there exists a range of linear charge density values, where our solutions may be 
stable with respect to small azimuthal perturbations. 

PACS numbers: 47.65.-|-a, 41.20.Cv, 47.20.Ma, 47.27.Wg 



I. INTRODUCTION 

Cylindrical jets are known to be unstable with respect 
to small surface perturbations because of the develop- 
ment of the Rayleigh instability caused by capillary ef- 
fects J. For electrically charged jets, electrostatic forces 
are an additional factor that determines the system be- 
havior. The Coulomb interaction of electric charges sup- 
presses the large-scale axial capillary instability. On the 
contrary, it can lead to the growth of nonaxisymmet- 
ric modes of disturbances which are stable for the un- 
charged jet (see and references therein) . In order 
to understand the main laws governing the behavior of 
a charged jet, it is important to define conditions when 
the mutual compensation of the electrostatic and cap- 
illary forces is possible, as well as the conditions when 
such compensation is impossible. Therefore, the neces- 
sity arises to determine the region of existence of stable 
solutions for the problem of the equilibrium configura- 
tions of the jet surface. 

In this paper we consider possible equilibrium shapes 
of a charged infinite jet of a conducting liquid (an elec- 
tric charge distributes itself over the liquid surface so that 
the electric field potential is constant everywhere inside 
of the conductor). In doing so, we restrict our analysis 
to the particular case of the azimuthal deformations of 
the initially circular jet (axial deformations can be sup- 
pressed by the longitudinal magnetic field). For this case, 
the instability is induced by the electrostatic forces, while 
the surface tension forces play a stabilizing role. 

We are presently only aware of a few nontrivial solu- 
tions of the classical problem in electrostatics, that is the 
problem of finding the stationary configurations of the 
charged surface of a conducting liquid (the flat surface, 
the circular cylinder surface, and the sphere belong to the 
trivial solutions). We should primarily mention the so- 
called Taylor cone. In Ref. j^, Taylor has demonstrated 
that the surface electrostatic pressure for a cone with an- 
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gle 98.6° is inversely proportional to the distance from 
its axis and, hence, can be counterbalanced by the capil- 
lary pressure (the sole exception is the cone apex where 
the force balance condition is violated). Recently, the 
first author of the present paper has shown |^ that for 
the case of plane symmetry, the problem of the steady- 
state shape of the free surface of a conducting liquid in 
an external electrical field is mathematically similar to 
the problem of the progressive capillary wave solved by 
Crapper . The analogy provides an easy way of solv- 
ing the electrostatics problem Q . Finally, the method of 
constructing exact solutions for the equilibrium config- 
urations of the charged two-dimensional drops was pro- 
posed in Ref. . It should be noted that apart from the 
transformations, this problem coincides with that of the 
stationary shape of a two-dimensional air bubble in a cir- 
culatory ambient flow. The solutions found by Zubarev 
for an arbitrary mode number n were independently 
given by Crowdy for n = 2 llOl and by Wegmann and 
Crowdy for n = 3,4,5... in the problem with the 
bubbles. The approaches developed in Refs. 0, 0, ^ 
turn out to be useful for the following analysis of possible 
configurations of charged jets. 

The article is made up as follows. In section II we give 
the equations defining the equilibrium configuration of 
the charged surface of a conducting liquid for the case of 
plane symmetry. It is shown that conformal transforma- 
tions allows us to reduce the investigation to the analysis 
of a nonlinear boundary- value problem on a half-plane for 
the Laplace equation. In section III, using the results of 
Refs. we obtain exact solutions for the jet con- 

figurations corresponding to the azimuthal mode num- 
bers n = 2,3,4... In section IV the equilibrium surfaces 
corresponding to our solutions are investigated. We for- 
mulate the conditions under which the surfaces become 
self-intersecting, and the jet splits into several separate 
jets. In section V we analyze the dependence of the jet 
surface deformation amplitude on the control parame- 
ter (linear electric charge density) for different azimuthal 
numbers. It turns out that we deal with a soft loss of sta- 
bility of the round jet surface for n = 2, 3, 4 (supercritical 
bifurcation) and with a hard loss of stability for n > 4 
(subcritical bifurcation). Section VI contains our conclu- 
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FIG. 1: The geometry of the azimuthaUy perturbed hquid jet 
is represented schematically. The jet cross-section is constant 
along the z-axis. 

sions and some remarks concerning conditions whereby 
the solutions obtained can play an important role in the 
jet behavior. 

II. INITIAL EQUATIONS 

Let us write the equations of electrostatics that de- 
scribe a stationary profile of the charged surface of the 
conducting liquid jet with constant cross-section along 
the direction of its motion (see Fig. The distribution 
of the electric field potential in the plane of the jet 
cross-section {x, y} is determined by the Laplace equa- 
tion: 

iPxx + Vyy = 0- 

It should be solved together with the condition that the 
conductor surface is equipotential, 

(^ = 0, 

and also the condition that the field of the charged con- 
ductor coincides at infinity with the field of a uniformly 
charged straight filament: 

tp ^ -Q\n{x^ + y'^), x'^+y^^oo, (1) 

where Q is the linear electric charge density of the con- 
ductor. 

The jet will be considered to move at the constant 
velocity v along the z-axis of the Cartesian coordinate 
system (a fluid is at rest in the system of coordinates 
moving with the jet). Then the equilibrium relief of the 
charged boundary of a conducting liquid is determined by 
the Laplace- Young stress condition, that is the balance 
condition for the electrostatic and capillary forces acting 
on the surface: 

(8^)-i(V^)2^o + «C' + P = 0, ^ = 0, (2) 

where a is the surface tension coefficient, and C is the lo- 
cal curvature of the surface. The constant p is expressed 



in terms of the jet velocity {v), hquid density (p), and 
the external (pe) and internal (pi) pressures: 

p^Pi-Pe- pv^ 12. 

It is apparent that an infinitely long cylindrical jet with 
a circular cross-section (i.e., a round jet) gives the trivial 
solution of the problem. In what follows, we will obtain 
its nontrivial solutions corresponding to the azimuthal 
deformations of the round jet surface. 

For convenience we convert to the dimensionless vari- 
ables: 

X ^ Q^{2'KaY^ X, y ^ Q'^{2TTa)~^ y, 
ip^2Qip, p^2na^Q~^p. 

By analogy with Rcf. we choose / = In jV^?] as a 
new unknown function, and the pair of conjugate har- 
monic functions ip and ip as new independent variables 
(the condition = const defines the electric field lines). 
The so-called complex potential w = ip — iip is an analytic 
function of the complex variable z = x + iy. The complex 
expression \n{—dw/dz) is also an analytic function, and, 
as a consequence, the real functions 

/ = Re ln(-du;/dz) = ln|V(^|, (3) 
9 = — Im ln{—dw/dz) = arctan{ipy/ipx) (4) 

are conjugate harmonic functions of the variables ip and 
-0. In particular, this implies that the function / satisfies 
the Laplace equation: 

fipip + fipTp = 0. (5) 

The boundary conditions for / can be derived from the 
expressions (QJ and We get 

U^pe-f + e^, (^ = 0, (6) 
f ^ LP, p^-oo, (7) 

where we have taken into account that the fluid surface 
curvature can be expressed in terms of the functions / 
and (p in the following way: C = — (/^ exp /)y=o. Since 
in the limit |2;| ^ oo we have w ^ —Inz for the complex 
potential and, consequently, a closed surface corresponds 
to changing ijj by 2tt, we add the condition for periodicity 
of / with respect to the variable ip, 

f{^,ij)^f{p,^ + 27T), (8) 

closing the system of equations for the function /. 

Thus, the problem of finding the steady-state pro- 
file of a cylindrical jet surface amounts to studying the 
boundary-value problem lO-® on the half-plane p < 0. 

III. EXACT SOLUTIONS 

A wide class of particular solutions of the equations 
©-(IHJl was obtained in Ref. 9]. They are given by the 
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formula 



, /l + a262e2"'^+2a6e"'^cos(nV'r 
In ( — rrr^r- ^ ■ „ , ^— f I + if, 



q2 _|. 52g2nv _ 2a6e"¥' cos(nV) 



where we put 



(9) 



In these expressions n is the azimuthal mode number 
(n = 2, 3, 4...) and I = ^/l — 4p is the parameter charac- 
terizing the amphtude of the surface deformation of the 
round jet. 

Let us construct the equihbrium surfaces correspond- 
ing to the solution © . It follows from the definitions © 
and Q that the inverse transformation from the vari- 
ables if and to X and y is determined by the relation 



exp(— / + iO)dw, 



(10) 



where is the angle of inclination of the electric field 
intensity vector to the abscissa direction. Taking into 
account that / and 9 are conjugate harmonic functions 
and thus the Cauchy-Riemann conditions are satisfied 
{Q^ji = fip and 9^ — —fip), we obtain from Eq. I^: 



f -i9 = In 



l-l 
2i 



2 In 



1 + a6e'"" 
a - fee""" 



Substituting this expression into (|10|l . we get: 

2i ffaf'-bV 2iaH(t" + b/a^) 



(l-l) 



-1- ab 



{I - 1) -I- ab) ' 



(11) 



where we have introduced the notation t = exp(— u;). 
The transformation corresponding to ifTTl) maps the unit 
circle |i| = 1 onto the free surface of a liquid. The alge- 
braic expressio n lllll re lates the results of Rcf. ^ and the 
results of Refs. \wL [ll| , where the conformal mapping to 
the exterior of the init disc was applied for the problem 
related to an analysis of the profile of a two-dimensional 
bubble in a circulatory ambient flow. 

Bearing in mind that f — dA, the boundary and, 
consequently, t = exp(i'!/)), we obtain from that the 
sought-for equilibrium surfaces are given by the following 
parametric expression: 



2ia^e''l'{e"''^+b/a^] 
{I - 1) (e™'^-|-a6) 



(12) 



where V plays the role of the parameter. The closed 
surface corresponds to the change in ip in the range < 
tp < 27r. 

Let us return to the real variables. Having separated 
the real part from the imaginary one in 112|) . we find: 



y = - 



a ^6cos {nij — ijj)+a'^b cos {nilj+ip) + {a'^ + b'^) costp 



{1 + a'^b^ + 2ab cos {n'ilj)){l- l)/2 



(13) 




FIG. 2: Typical (superposed) cross sections of the charged 
jet of a conducting liquid for n = 2 and parameter values 
I = 1.86, 1.94, 1.88, 2. The cross section areas are normalized 
to a constant. 



a ^bsin {ml> — 'il>) — a^b sin (n'ljj+'ijj) — (a^ + 6^) simp 
^ " (1 + a2&2 + 2abcos{niP)){l - l)/2 ' 

(14) 

The formulas (|13|l and H14|) represent a family of exact 
two-parametric solutions for the equilibrium shape of a 
charged jet of a conducting liquid. To the best of our 
knowledge, these solutions for the jet configurations have 
not been considered so far. 

The particular solutions (|13|l and H14(l of the initial 
equations relate to the case p 0- Note that for p = it 
is possible to find the general solution of Eqs. ISj-JSJ. It 
has the following form: 

/ = In (l - c^) - In (c^ - 2ce-'^cos ip + e"^'^) - (p, (15) 

where the constant c satisfies the inequality < c < 1. 
This solution is 27r-periodical with respect to -0, so that 
it corresponds to the azimuthal number n = 1. The 
equilibrium surface is given by the parametric expression 



1 



1 -c2 



lip 



2ciTp) . 



One can readily see that this surface has the self- 
intersections for arbitrary c. Thus, the solution l|15|) with 
the mode number n = 1 is physically meaningless. 



IV. CONDITIONS OF THE JET SPLITTING 

The solutions (|13|l and (|14|l obtained in the previous 
section enable us to find exact critical values of the linear 
charge densities required (i) for the onset of the azimuthal 
instability of the round jet and (ii) for the jet splitting. 

In the limit / — > n expressions (|13|l and H14|l for the 
equilibrium shape of a charged jet of a conducting liquid 
define circles of radius 2/{n + 1), which correspond to 
the unperturbed state of the jet, namely, to the round 
jet. With decreasing the parameter I, the jet surface is 
deformed. At certain n-dependent critical values In of 
the parameter Z, the region occupied by the liquid ceases 
to be simply connected, and the jet splits (see Fig. 
For 1 < I < In the solutions are physically meaningless 
so that for fixed n the set of the problem solutions cor- 
responds to the interval In < I < n. 
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FIG. 3: The cross sections of the jet for n — 3 and parameter 
values I = 2.53, 2.59, 2.78, 3. 




FIG. 4: The cross sections of the jet for n — 4 and parameter 
values I = 3.19,3.35,3.8,4. 

We now determine the critical values of the parameter 
I. For n — 2 the condition of the surface self-intersection 
has the form 

a; = 0, V = 71-/2, 

that corresponds to the singularity at the point x = y = 
0. It follows from this condition that I2 is a root of the 
quadratic equation 

7^^-6^2-13 = 0. 

Only one solution of the equation, I2 = 13/7 « 1.86, 
meets the requirement 1 < I2 < 2. The jet splits into two 
equal parts at this value of the parameter I (see Fig. ??). 

For n > 2 the condition of the self-intersection of the 
surface can be written as 

a; = 0, dx/dip = 0- 



Hence it follows that the critical values oil, i.e., the quan- 
tities In, are determined by the set of equations: 

a{n + 1) (a^ + b^) cosi/) 
= 2h{n — 1) [n sin(n-i/;) sin ijj + cos(nV') cos -0] , (16) 

a{n -t- 1)^ (a^ -1- h^) smi) 
= 2b [2nsin(n-0) cos V' — {n^ + 1) cos{mp) sin-0] . (17) 

It is easy to see that, dividing (|17|l by H16|l . we can elim- 
inate the parameter I from these equations. As a result 
we obtain the equation for -0 at the point of the curve 
self- intersection: 

tan{nip) [(n^ - l)tan^?A ^ 2] + 2n tanip = 0. 

With the help of de Moivre's formula, this trigonometri- 
cal equation can be brought to the algebraic form: 

[{n^-l)h^-2] lm{l + ihy'+ 2nhRe{l+ih)" ^ 0, (18) 

where we put h — tan?/;. Solving this equation of the nth 
order, we can define the value of h and, consequently, of 
■0 for any n. At given 0, the sought-for critical values of 
the parameter I can be determined from the relation H16|l . 
which can be transformed into the quadratic equation 
with respect to the unknown quantity /„. 

For n — S the expression p8|l transforms into the triv- 
ial equation /i^ = 1, from which it follows that the self- 
intersection takes place at ip = it / 4. For this value of the 
parameter -0, the condition (|15|l reduces to solving the 
equation 

17/^-18/3-63 = 0. 

Its root satisfying the condition 1 < /a < 3 is /s = (9 -I- 
24-\/2)/17 ~ 2.52. If / = /s, the volume occupied by liquid 
loses its simple connectivity and one jet splits into four 
unequal parts. This situation is illustrated in Fig. |21 

For n = 4 (see Fig.EJ, we obtain from = 5/13, 

and H16|l transforms into the equation: 

443 486/4 - 2957 = 0. 

It follows herefrom that U = (243 + 370VT0)/443 w 3.19. 
For the next azimuthal number, n = 5, we deduce from 
Eq. IIHI): 

3/1"* - 24/1^ 5 = 0, 

and hence = A± ■y/43/3. Substituting this value of 
the parameter h into H16|l . we get: /s w 3.85. In a similar 
manner, we can find that /g ~ 4.55 and /s ~ 5.85. 

Thus, we have defined the values of the parameter / 
for which the solutions of the problem of the equilibrium 
configurations of the jet surface exist. Let us now con- 
sider a jet with given characteristics (the cross section 
area S and the surface tension coefficient a) and deter- 
mine the linear charge densities Q, which correspond to 
the allowable values of /, i.e., to the interval /„ < / < n. 
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For the family of the solutions (|13|l and H14() , the area S 
can be easily found with the help of the Green's formula 
(compare with Ref. |0|): 



1) 



An) 



2iT/n 


where the dependence (|12f) of the complex variable z on 
have been used. Returning to initial dimensional quan- 
tities, we get 



S 



{{I + 1)^ - An) 
7ra2(Z2-i)2 ■ 



Solving this relation with respect to Q, we arrive at the 
dependence of the linear charge density on the surface 
tension a, the jet cross section area 5, and the steady- 
state solution parameters / and n: 



Q = 



-,1/4 



An 



(19) 



It enables us to determine the critical values of the charge 
density. 

In the problem under consideration, we can define 
two critical charge density values for every azimuthal 
wavenumber n. The first critical density value Qn cor- 
responds to the threshold of linear instability of the jet 
with round cross section. It can be calculated from the 
formula (|19|l . where we must take / = n (recall that for 
I — n the expressions ()13f) and H14|l define circles): 



Qn = [An + l)^a^S] 



1/4 



The critical charge density is seen to grow monotonically 
with n. Its minimum value corresponds to n = 2: 

Q2 = (9™2 5)'/'«2.31ai/25i/4. 

If the linear charge density Q of the jet exceeds this value, 
an electrohydrodynamic instability of the jet surface will 
develop. At the initial stage of the process, this leads to 
the elliptic deformation of the jet cross-section. The dis- 
turbances corresponding to the modes with lesser spatial 
scales, n > 2, can grow only at larger values of the linear 
charge density. 

The second critical density value Qn corresponds to 
the situation when the volume occupied by the jet loses 
its simple connectivity and the jet splits into two or more 
separate jets (see Fig. We can find it substituting 

I = In into ifT^ll : 



{In 



1 



An 



1/4 



The charge density takes the minimal value for n — 2, 
whence it follows that the jet splitting into two approx- 
imately equal parts can be considered as the most prob- 
able scenario of the jet disintegration. Since I2 = 13/7, 



the exact minimal value of the second critical density is 
given by the following expression: 

Q2 = (l800 7ra25/49)^^« 3.28 ai/25.1/4 

Note that the critical densities Q„ can be found from 
the linear analysis of the stability of the round jet surface 
(see, for example, Ref. 0, 0)- In so doing it is not 
necessary to know the exact solutions for the stationary 
jet shape. However we must know them for determining 
the conditions of the jet splitting or, what is the same, 
for finding the second critical charge densities. 



V. STABILITY ANALYSIS 

Now consider the dependence of the jet surface defor- 
mation amplitude on the linear electric charge density. 
Such an analysis will allow us to draw some qualitative 
conclusions concerning the stability of the solutions ob- 
tained. 

It is convenient to take Q2 as a unit of linear elec- 
tric cha rge d ensity and the radius of the unperturbed jet 
i?o — •\/5'/7r as a unit of length. This implies introducing 
both dimensionless charge density q and the deformation 
amplitude of the jet surface r: 



-,1/4 



r{nj) 



Rn 



Ro 



9((; + l)2 -4n) 



2V^^2— p + (1-1)^:^^2—1 
v/(n2 - 1) ((/ + 1)2 - An) 



- 1, 



where i?niax = y\ii>=o is the maximum distance between 
the jet axis and its surface. These relations give the de- 
pendence of r on (J for different n in the parametric form 
{I plays the role of the parameter). The relevant plots are 
presented in Fig. [S] The straight line r = in the figure 
corresponds to the unperturbed state of the system, i.e., 
to the round jet. 

One can notice that the deformation amplitude mono- 
tonically increases with the charge density for n ^ 2,3, A. 
This situation corresponds to the soft loss of stability 
of the round jet. Indeed, it is clear that if q < Q/Qn, 
then the round jet is stable with respect to small sur- 
face perturbations with azimuthal wavenumbers greater 
than or equal to n. Bifurcations occur when the condition 
q = Q/Qn holds (supercritical bifurcation for n = 2, 3, 4). 
It is seen from the figure that the side branches corre- 
sponding to our exact solutions fork from the strai ght line 
r = 0. It should be noted that the authors of Ref. [Tj, in 
terms of the present paper, have plotted the dependence 
of g2 on p that also indicates the manner in which the 
solution branches bifurcate. 

As the trivial solution r = is unstable for (7 > 1, the 
free energy of the system can have a minimum on these 
branches only. This suggests that the surface modes with 
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FIG. 5: The dimensionless amplitude of the jet surface defor- 
mation r as a function of the normalized linear charge density 
q for different azimuthal numbers n. The straight line r = 
corresponds to the unperturbed state of the system, i.e., to 
the round jet. 

small azimuthal wavenumbers are excited in a soft regime 
and, at least for a small overcriticality, our exact solutions 
can be stable with respect to small pertm'bations that do 
not violate the problem symmetry. 

For n > 4, the surface deformation amplitude r de- 
creases with increase in q in the vicinity of the branch 
points (see Fig.jSJ. Such a dependence of the amplitude 
r on the control parameter q corresponds to the subcriti- 
cal bifurcation. Then the potential energy of the system, 
i.e., the sum of the surface and electric field energies, has 
a maximum on the side branches. It immediately fol- 
lows that our solutions having azimuthal mode numbers 
n > 4 are unstable with respect to small changes in the 
amplitude r. 

Since both the first and second critical charge densities 
are minimal for the large-scale azimuthal mode n = 2, it 
is the surface mode that will define the jet behavior. If 
the linear charge density Q exceeds the value Q2, then 
the surface of the cylindrical jet with a round cross sec- 
tion becomes unstable. As the instability regime is soft, 
then, for a small subcriticality, a new stable state cor- 
responding to our stationary solutions with n = 2 ap- 
pears, and the round jet transforms to an elliptic one. 
Unfortunately, the analysis of the balance conditions for 
the forces acting on the jet surface cannot provide the 
answer to the question as to whether our solutions are 



stable for all permissible values of the electric charge 
density, Q2 < Q < Q2 (we discuss the stability with 
respect to surface perturbations that do not violate the 
problem symmetry) . In all probability, they are unstable 
for a sufficiently large Q. The point is that Q3 < Q2 
and Q4 < Q2, and, consequently, the functional of the 
jet potential energy can have several extremums corre- 
sponding to different n at given Q. Moreover, as evident 
from Fig. 12 the jet cannot be stable if the charge density 
value is close to Q2- An arbitrary small deformation of 
the surface can lead to the rupture of the neck. Thus the 
condition Q > Q2 may be considered as the sufficient 
condition for splitting the jet into two. 



VI. CONCLUDING REMARKS 

In the present work we have obtained the two- 
parameter family of the exact solutions of the classical 
problem in electrostatics, namely, the problem of finding 
the equilibrium configuration of a charged jet of a con- 
ducting liquid. The approach applied is based on the con- 
formal mapping of the region outside the jet to the half- 
plane, that have restricted our consideration to the case 
of the plane symmetry of the problem, when all quanti- 
ties depend only on two variables x and y (see Fig. 
Because of this, all the solutions obtained correspond to 
the azimuthal deformations of the jet surface, whereas 
the deformations in the direction of the jet axis have not 
been considered. At the same time, it is just the lon- 
gitudinal instabilities (the varicose and sinuous modes) 
that determine the behavior of a charged jet in the gen- 
eral case 0. Nevertheless, if the electrohydrodynamic 
instability of the liquid cylinder is suppressed in the di- 
rection of the z axis, our solutions can play a dominant 
role in the jet behavior. For instance, the growth of the 
axial disturbances can be stabilized by the magnetic field 
directed along the jet axis. It is known 0] that the tan- 
gential magnetic field retards the development of the sur- 
face instabilities which bend the field lines (in particular, 
this phenomenon is used to confine the plasma). If the 
magnetic field is sufficiently strong, the only azimuthal 
instability of the jet surface will develop. As it was dis- 
cussed in the previous two sections, such instability leads 
to the jet splitting. 

This study was supported by the President of Rus- 
sian Federation (Project No. MK-2149.20042), the Foun- 
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